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MODELING LUNAR VOLCANIC ERUPTIONS. R. M. Housley, Rockwell Inter- 
national Science Center, Thousand Oaks, California 91360. 

From the time the first returned lunar samples were examined, it has been 
known that early in its history the moon underwent a period of volcanism, 
during which a large quantity of volatile poor basaltic magma was erupted. 

The nature of the volcanism and the influence of outgassing and volatile 
transport in the lunar vacuum have been extensively debated ever since. 

Although photogeologists interpreted features on several areas of the 
moon as pyroclastic deposits, it was not until after the Apollo 17 mission 
that pyroclastic materials were convincingly identified in the returned lunar 
samples. The case for so identifying the Apollo 17 orange and black glasses 
and the Apollo 15 green glass has since been systematically argued by Heiken 
et al. (1) on the basis of chemical, petrographic and photogeo logic evidence. 

A plausible explanation for the puzzling fact that these glasses do not 
resemble terrestrial pyroclastic materials has been offered by O'Keefe (2). 

He pointed out that the frothy nature of terrestrial pyroclastics may depend 
on the fact that the surface tension of the melt increases as water is 
vaporized from it. 

A plausible though modest source of volatiles that is compatible with 
the geochemical constraints imposed by the reduced nature of the lunar basalts 
has recently been suggested by Sato (3). He showed that all available data 
can be interpreted by assuming chat about 100 ppm of finely divided graphite 
was present in the magma at depth. Above about 3 km, this graphite would 
react with the melt to form bubbles containing CO and CO 2 and leave finely 
divided grains of Fe metal in the melt. This avoids any problem associated 
with the nucleation of bubbles in a volatile poor magma and could help explain 
the different physical characteristics of lunar and terrestrial pyroclascics . 

Thus, the stage is set for serious efforts to model lunar volcanic 
eruptions. We proceed in three steps, considering first, the transport of 
magma to the lunar surface; then the nature of the transient lunar atmosphere 
produced by the gas released; and finally, the processes that could occur in 
the vicinity of the point of eruption. 

As a starting point for modeling magma transport, we assume that the 
volcanic plumbing is similar for lunar volcanos and for terrestrial oceanic 
volcanos. From a visual inspection of recent Hawaiian craters, by allowing 
for near surface flaring, we obtain ar order of magnitude estimate of about 
10 ra for the typical diameter of terrestrial magma conduits. 

It is important to note that in a conduit of this size any reasonable 
estimates of flow velocity and magma viscosity yield very high Reynolds 
numbers, implying that the flow must be highly turbulent. In this case, 
the average flow velocity depends strongly on the roughness of the conduit, 
but only weakly on the viscosity (4). Thus, we reach the surprizing conclu- 
sion that the much lower viscosities of lunar magmas in comparison to their 
terrestrial counterparts will not lead to substantially increased eruption 
velocities . 

In both lunar and terrestrial volcanos, the ultimate driving force 
responsible for magma migration is the difference in average density between 
the fluid and the rock through which it moves. Although the near surface 
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behavior of terrestrial magmas is complicated by the large volatile component 
present, we may roughly check whether our estimated conduit diameter is realis- 
tic by comparing calculated eruption characteristics with the recorded behavior 
of Hawaiian volcanos (5). Assuming the maximum reasonable density contrast 
of 10%, a viscosity of 145 poise (1), and an average surface roughness of 0.5 m 
in a 10 m diameter conduit, we can use the hydraulic engineering data graphi- 
cally presented in (4) to calculate an average flow velocity of 16.5 m/sec, 
an eruption rate of 4.7 x 10^ m 2 /hr, and a Reynolds number of 3.4 x 10^. 

The eruption rate of 4.7 x 10° m 2 /hr is about a factor of 5 high for 
typical recent Hawaiian eruptions (5), probably because the assumed effective 
pressure head is too large. 

With the same conduit diameter and density contrast and a viscosity of 
10 poise (1), the calculated lunar eruption rate is 1.9 x 10° m J /hr, with a 
Reynolds number of 2 x 10 5 and a flow velocity of 6.8 m/sec . 

From he Fe metal content of lunar basalts, we can calculate the amount 
of CO they must have evolved (3). Using 0.15 wt% Fe metal as determined by 
Morris and Gose (6) for 74001 glass, this is 0.075 wt%. Thus, our hypothetical 
eruption brings CO to the lunar surface at the rate of about 1.2 x 10^ gm/sec. 

As a first approximation, we assume that at modest distances from the 
source this gas expands hemispherically . A limiting velocity of about 1.5 km/ 
sec will be reached when all the thermal energy is converted into kinetic 
energy of flow. The density p and velocity v at any radius r and mass flow 
rate m are related by 

2:rr 2 pv = m 

At distances where v is close to its limiting value with the above mass 
flow rate 


P = 


1.3 x 10 -4 


R 


gm 


cm 


where R is the radius in m. The °C pressure corresponding to this p is 


„ 1.04 x 10 _L ^ „ 7.9 

P = r atoms, or P = 


10 


R 2 


— corr 


_o 

Thus, we see that the pressure only falls off to 10 torr at about 100 km 
from the eruption site. Yet this formula suggests that pressure nay remain 
quite modest near the eruption site, although a more complete calculation is 
needed to describe this region. 

Fountaining during terrestrial volcanic eruptions appears to be caused by 
the evolution of gas from the magma in the conduit as it moves toward the 
surface. The observed fountaining, in fact, requires the evolution of only 
a fraction of the tota. gas present. 

Our estimated molar abundance of CO in lunar magmas is about a factor of 
10 lower than the volatile content of terrestrial magmas. However, it may 
more completely evolve in the conduit and on an equal molar basis will occupy 
6 times the volume because of the lower pressure in the lunar conduit. Taking 
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into account these facts, we have estimated that lunar magmas could reasonably 
fountain to a height of about 1 km. 

At the point of eruption the gas would occupy considerably more than half 
the volume and would still be under sufficient pressure to expand violently 
separating much of the material into small droplets with velocities of the 
order of tens of m/sec. 

This modeling effort clearly can be refined and extended in many ways. 

For example, the ability of the expanding gas cloud to transport solid particles 
and perhaps to generate lunar dust storms can be investigated. 

At the present stage, it suggests several important experiments. The 
kinetics of graphite reaction with Fe silicate melts should be investigated. 

The partitioning of trace volatiles between gas and liquid should be studied. 

And the behavior of a dry vesicle ladden melt on sudden decompression should 
be checked. 

It also strongly suggests the following conclusions: 1) The difference 

in viscosity between lunar and terrestrial magmas is of minor importance in 
comparing their eruption characteristics. 2) Magma may fountain to heights 
of the order of 1 km during lunar eruptions. 3) Violent expansion of CO 
bubbles at the instant of surface arrival can lead to the formation of 
abundant ~ 100 um glass balls and their distribution over a radius of hundreds 
of meters. 4) A major fraction of the more volatile elements in the magma 
may be in the bubbles and, hence, escape into the transient atmosphere. 

5) The expansion of the transient atmosphere will distribute volatile elements 
essentially over the entire moon. 

(1) Heiken G.H., McKay D.S., and Brown R.W. (1974) Geochim. Cosmochim. Acta 38., 

p. 1703-1718. 

(2) O'Keefe J.A. (1976) Science 194 , p. 190. 

(3) Sato M. (1976) Proc. Lunar Sci. Conf. 7t'n, p. 1323-1344. 

(4) Standard Handbook for Mechanical Engineers, Seventh Edition, Baumeister T., 

Editor, McGraw-Hill, Section ?• 48-62. 

(5) MacDonald G.A. and Abbott A.T. (1970) Volcanos in the Sea, U. Hawaii Press. 

(6) Morris R.V. and Gose W.A. (1977) Proc. Lunar Sci. Conf. 3th, p. 3113-3122. 
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FLAMELESS ATOMIC ABSORPTION STUDIES OF VOLATILE TRACE METALS IN 
APOLLO 17 SAMPLES. E. H. Cirlin and R. M. Housley, Rockwell International 
Science Center, Thousand Oaks, California 91360. 

We report further developments in our study of vapor transport processes 
using a flameless atomic absorption technique (1,2) to analyze the volatile 
trace metals Pb and Cd in orange and black glass core samples 74002/ 74001 and 
nearby surface samples 74220 and 74241. 

Many workers have show, that volatile elements are enriched in the orange 
and black glasses and that, in general, the concentrations of these volatile 
elements are surface correlated (3,4). Previous workers have provided strong 
arguments supporting a volcanic origin for these glasses and have proposed 
lava fountain mechanisms (5,6). 

Our Pb results on dry sieved core tube sample 74001-105 from the depth 
of 39-39.5 cm also shows that the concentration of Pb increases as the grain 
size decreases (Fig. 1). A scanning electron microscope study (Fig. 2) of 
this black glass shows that most of the small grain size material is made up 
of broken pieces. This provides strong evidence that the volatile trace 
metals were deposited after the glass was broken. The glass balls could 
have been broken either upon landing on the lunar surface during the initial 
deposition from the volcanic eruption or by shock during the recent Shorty 
Crater event. 

Reed and his coworkers (7) and Heiken and McKay (5) both proposed 
separate lava fountaining mechanisms to account for the formation of the 
orange and black glasses. 

The total Cd (Fig. 3), total Pb and PbCTo (fig* A) concentration in 
core tube samples 74002/74001 decrease as core tube depth increases. Our 
study also shows that a considerable amount of Pb in the surface sample 
74220 and in the samples near the surface 74001-1069 is in the form of PbCl.2* 

The presence of orange glass down to a depth of about 25 cm in the 
upper core section 74002 and black glasses for the rest of 74002 and all 
of 74001 was noticed by Schmitt (8). 

If these orange and black glasses were formed during different lava 
fountains as suggested by Reed ei at. (7) and Heiken et ai. (5), then dis- 
continuity of the concentration of volatile elements at about 25 cm depth 
would be expected. Our data shows a continuity in concentration of volatile 
metals throughout the core. 

Our data suggest that either the volatile trace metals were completely 
redistributed by the Shorty Crater impact or that the original strata, if 
formed by a lava fountain, and assuming a single fountain mechanism, was 
totally preserved. Since it seems very unlikely that the original strata 
are well preserved at' the rim of Shorty Crater where the samples were 
collected, our results seems to suggest that volatile trace metals were 
redistributed and that PbC^* being volatile, was transported to the upper 
surface as vapor from below while the lower strata slowly cooled. 

The Pb concentrations in samples 74220 (2.29 ppm), 74241 (2.21 ppm) 
and the near surface sample 74001-1069 (2.40 ppm) are about the same, 
indicating a single exposure to Pb. 
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(1) Cirlin E.H. ami Housley R.M. (1977) Lunar Science VIII, p. 184-186. (2) Cirlin E.H. and 

llousloy R.M. (1977) Proc. Lunar Sci. Conf. 8th, p. 3931-3940. (3) Krahenbuhl U. Grutter A., 

Von Gunten H.R., Meyer G., Wegmuller F. , and Wyttenbach A. (1977) Proc . Lunar Sci. Conf. 8th, 
p. 3901-3916. (4) Wasson J.T., Boynton W.V., Kallemeyn G.W. , Sundberg L.L. and Wai C.M. (1976) 

Proc. Lunar Sci. Conf. 7th, p. 1583-1595. (5) llelkcn G. and McKay D.S. (1977) Proc. Lunar Sci. 
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AN XPS STUDY OF APOLLO 17 ORANGE AND BLACK GLASSES. R. W. Crane 
and R. M. Housley, Rockwell International Science Center, Thousand Oaks, 
California 91360. 

We report results of a detailed study of the surface chemistry of grain 
size fractions dry sieved from Apollo 17 black glass sample 74001,105. This 
sample from a depth of 39.0-39.5 cm in the double drive tube is free from 
agglutinitic contamination. T.n its finer fractions broken surfaces strongly 
predominate over original spherical surfaces. 

We have studied a 75-1000 ym sample cemented to a ~ecessed sample holder 
with colloidal graphite and two < 20 ym samples, one prepared according to 
our normal procedure by sedimentation from ethanol, and one packed dry into 
a recessed sample holder. These latter two samples gave indistinguishable 
results indicating that our data are not influenced by the method of sample 
preparation. 

Both size fractions showed high surface concentrations of Zn, Na, and S. 

No other minor or trace elements were present in comparable amounts. 

Although the 75-1000 ym fraction presents more than 50% unbroken surfaces, 
it only showed concentrations about a factor of two higher than those in the 
< 20 yn fraction, which presents at most a few % unbroken surfaces. Therefore, 
considerable redistribution of volatiles must have taken place after the grains 
were broken. This could have been either in the fountaining process itself or 
as the result of a subsequent impact. 

In Figs. 1 and 2 we present raw peak intensity versus sputtering time 
profiles for Zn and Na in the two size fractions. The s, uttering rate for 
a flat Si02 surface under the conditions used is 0.7A per minute. It appears 
that although the surface enrichments are chin in both cases that they are 
somewhat thicker on the unbroken surfaces. 

When we initially ran the samples after pump down at room temperature, 
about half the S was in an oxidized form. Counting both the oxidized and 
sulfide forms, there still was about a factor of two too little S present 
to correspond with the observed Zn. Before sputtering, we heated the samples 
to 300°C for an hour. This removed the majority of the oxidized S in agree- 
ment with mass spectroscopic results previously reported by Simoneic at. (1). 

We feel the data suggest that the Zn was originally transported and 
deposited as ZnS and that much of the S has been oxidized and lost since the 
samples were returned. 

(1) Simoneit B.R., Neil J.M., Wszolek ?.C., and 3urlingar.e A.L. (1976) In 

Lunar Science VII, p. 815-816. The Lunar Science Institute, Houston. 
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a of Zn and Na on 75-1000 pm black glass 74001 as a function 
tial rise is due to removal of carbon contamination. 
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TEMPERATURE DEPENDENCE OF THE QUANTITATIVE FERROMAGNETIC RESONANCE 
OF SIZE SEPARATED FINES SAMPLES. I. B. Goldberg, R. M. Housley, and H. R. 
Crowe, Rockwell International Science Center, Thousand Oaks, CA 91360 

The excess iron metal in lunar fines samples has been shown to be re- 
lated to their surface exposure ages (1) . Recent techniques permit the iron 
to be quantitatively determined (2), and these aethods are currently in use to 
measure the relative surface exposure ages of core samples (3-5). 

The line shapes and widths of the FIR absorption vary among lunar fines 
samples for reasons which are not understood. Linewidths and shapes as well 
as specific intensities also differ in different size fractions from the sane 
sample. Puzzling temperature and frequency dependances of the linewidths and 
intensities have also been reported in the past. We are studying the tem- 
perature and frequency dependance of the FRM absorption in size fractions of 
fines samples in the hope that we can understand the mechanisms that control 
shape and intensity and can in turn use this information to expand our under- 
standing of the evolution of the samples. 

The temperature dependencies of the peak-to-peak linewidths of size 
separates of sample 12001 are shown in Fig. 1. At room temperature, the line- 
widths are equivalent, but as the temperature is lowered, the linewidths of 
the smaller grains increase much more rapidly than that of the larger grains. 
The derivative and computed absorption curves shown in Fig. 2 show that for 
the l-2pm separate, this linewidth increase is not merely due to the change 
in the anisotropy field. While this increase in the linewidth may be due to 
an increasing contribution from a broader resonance component, the maxima of 
the linewidth at 20-30K for the 2-5um and 20-30um size separates parallel 
valence-exchange or impurity broadening of FMR spectra (6). This broadening 
was shown to occur due to surface oxidation of Ni-Fe films (7). 

Double integral values of the l-2pm size fraction increase by 39" bet- 
ween 270K and 10R, while those of the 2-5pm size range increase 2S”, and those 
of the 20-30ym size fraction increases by 21". These results indicate that 
the smaller grains contain a greater proportion of spheroidal iron which is 
superparamagnetic at high temperatures . We observed only a monotonic increase 
in the double integral values, and no decrease below 200K 3S reported by 
Griscom et. al. (8). 

Since all iron particles are saturated above about 21K0e, and spherical 
particles are saturated above about 7K0e, quantitative 35 GHz determinations 
were carried out to estimate the content of mult i-denain Fe® particles rela- 
tive to the 9.5 GHz determinations. The results of these measurements are 
shown in Table I. 
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TABLE 1 

QUANTITATIVE FMR MEASUREMENTS AT 


9.5 AND 34 

GHz FOR 12001 

SAMPLES (294K) 


Size Range Weight 

(pm) (rag) 

7. Fe° at 
9.535GHz 

Z Fe* at 
33.982GHz 

Relative 

Increase 

0.5-2 0.079 

1.3-4 

1.72 

1.28 

10-2 0.099 

0.553 

0.751 

1.35 

35-45 0.124 

0.350 

0.525 

1.50 


All samples showed an apparent increase in the amount of Fe° when 
measured at 35GHz. The finer size fractions exhibited a smaller increase 
indicating that they contain smaller spheroidal iron particles. 

During *n- e measurements it became ap_ .rent that as the amount of 
sample used incr^ared, the percentage of Fe° appeared to decrease. This 

was attributed to loading of the microwave cavity, as shown in Table II. 


TABLE II 


Sample 

wt 

35CH 

Z 

9 

. 5GHz 

( 45 m) 

(mg) 

AHpp 

Tfe T 

AHpp 

%Fe J 

6570i,13 

0.063 

672 

• 17 

570 

0.825 


0.552 

748 

87 

564 

0.471 


3.00 

- 

- 

556 

0.470 

67701,26 

0-118 

687 

0.195 

554 

0.105 


0 951 

725 

0.140 

547 

0.115 


2.20 



564 

0.124 

66031 

0.093 

694 

0.373 

578 

0.323 


2.25 



572 

0.390 


This anomalous decrease in the X a band spectra results from loading the 
microwave cavity so that the intensity is no longer linear. Several results 
report large variations of the linewidth of the 35GHz resonance while that of 
the 95GHz resonance is constant. Using samples smaller than about 150ug at 
35GKz, (depending on the amount of Fe°), the linewidths are constant to within 
215 Oe and the percentage of Fe° is greater than that measured at 95GHz. 

REFERENCES 

(1) Housley, R.M., Grant, R.W. , and Paton, N.E. (1973) Pro'*. Lunar Sci. Conf. 
4th, p. 2737-2749. 

(2) Housley, R.M. , Cirlin, E.H. , Goldberg, I.B., Crowe, H.R., Weeks, R.A., 
Perhac, R. (1975) Proc. Lunar Sci. Conf. 6th, p. 3173-3136. 

(3) Cose, W.A. and Morris, R.V. (1977) Proc. Lunar Sci. Conf. 8th, p. 2909- 
2928. 
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